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Abstract. We propose an approach to the theory of Riesz transforms in a framework 
emerging from certain reflection symmetries in Euclidean spaces. Relying on Rosler's 
■ construction of multivariable generalized Hermite polynomials and Hermite functions 

O . associated with a finite reflection group on M. d , we define and investigate a system of 

Riesz transforms related to the Dunkl harmonic oscillator. In the case isomorphic with 
the group Z2 it is proved that the Riesz transforms are Calderon-Zygmund operators in 
the sense of the associated space of homogeneous type, thus their mapping properties 
follow from the general theory. 



< 

u 

1. Introduction 

In [6] the authors developed a unified approach to the theory of Riesz transforms in 
the setting of mult i- dimensional orthogonal expansions with a rather general second order 
differential operator as the underlying "Laplacian". It is remarkable that this approach 
also works within the framework of differential- difference operators on M. d related to finite 
reflection groups, the setting intimately connected with the Dunkl theory, which has gained 
a considerable interest in various fields of mathematics as well as in theoretical physics 
qq ■ during the last years. 

Given a finite reflection group G C 0(M. d ) and a G-invariant nonnegative multiplicity 
function k : R — > [0, 00) on a root system R C M d associated with the reflections of G, the 
Dunkl differential-difference operators T!f, j — 1, . . . , d, are defined by 

c3 ■ 

T*f(x) = djfix) + Yl / e C\R d ); 

/3eR+ \Pi x l 

here dj is the jth partial derivative, (•, •) denotes the Euclidean inner product in M d , 
R + is a fixed positive subsystem of R, and ap denotes the reflection in the hyperplane 
orthogonal to (3. The Dunkl operators Tj, j = 1, . . . , d, form a commuting system (this 
is an important feature of the system, cf. [2]) of the first order differential-difference 
operators, and reduce to dj, j = 1, . . . , d, when k = 0. Moreover, are homogeneous 
of degree —1 on V, the space of all polynomials on M. d . This means that T*V m C V m -\, 
where meN = {0,l,...} and V m denotes the subspace of V consisting of polynomials of 
total degree m (by convention, V-\ consists only of the null function). 



Key words and phrases: Dunkl operators, Riesz transforms, Calderon-Zygmund operators. 
Research of both authors supported by MNiSW Grant N201 054 32/4285. 

1 



2 



A. NOWAK AND K. STEMPAK 



In Dunkl's theory the operator 

d 

plays the role of the Euclidean Laplacian (in fact A comes into play when k = 0). It is 
homogeneous of degree — 2 on V and symmetric in L 2 (M. d ,Wk), where 

if considered initially on C^° G (M. d ), the space of C°° functions on IR d with compact support 
outside the union of all the hyperplanes Hp, (3 G R + ; here Hp denotes the hyperplane 
passing through the origin and orthogonal to (3. Note that is G- invariant. 

In this article we propose a definition of Riesz transforms associated to the operator 

Lk = — Afc + ||a;|| 2 , 

which, due to the harmonic confinement \\x\\ 2 , we call the Dunkl harmonic oscillator 
(it becomes the classic harmonic oscillator —A + ||x|| 2 when k = 0). In the case of 
a reflection group isomorphic to the group we study L p mapping properties of the 
introduced Riesz transforms in detail. In the general case it occurs that (or rather its 
self-adjoint extension has a discrete spectrum and the corresponding eigenfunctions 
are the generalized Hermite functions defined and investigated by Rosier [10] . Then, under 
some mild assumptions, the formal definition Rj = 5j{Ck)~ 1 ^ , rewritten in terms of the 
related expansions, delivers an L 2 -bounded operator. 

The results of the present paper naturally extend those established in [13] by J.L. 
Torrea and the second author. On the other hand, the results are closely related to the 
authors articles [7], [8], where Riesz transforms for multi-dimensional Laguerre function 
expansions were defined and thoroughly studied. For basic facts concerning Dunkl's theory 
we refer the reader to the excellent survey article by Rosier [10]. There, one can also find a 
discussion (see [TUl Section 3]) and extensive references concerning applications of Dunkl's 
theory in mathematical physics. 

The organization of the paper is the following. In Section [2] we define Riesz transforms 
in the context of the harmonic oscillator based on a general Dunkl operator. Section [3] 
introduces the particular Dunkl setting related to the group 7L\. Here, apart from gath- 
ering basic facts, we establish a new result (Theorem 13.11) concerning the heat semigroup 
maximal operator. In Section H] the 7L\ Riesz-Dunkl transforms and the relevant kernels 
are defined, and the main results of the paper are stated (Theorem 14.31) . As a typical ap- 
plication of L p -boundedness of Riesz transforms, a priori L p -bounds in the context are 
then derived (Proposition 14.41) . Finally, Section [5] is devoted to the proofs of all necessary 
kernel estimates related to the setting. 
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Throughout the paper we use a fairly standard notation. Given a multi-index n G N d , 
we write \n\ — n\ + . . . + n^, n\ = n^. ■ . . . ■ and, for x, y G M. d , x n = x™ 1 ■ . . . • x^ d , xy = 
(xiyi, . . . ,Xd,yd)'i \\ x \\ denotes the Euclidean norm of x G M d , and tj the jth coordinate 
vector in M. d . Given x G M. d and r > 0, B(x,r) is the Euclidean ball in M. d centered at x 
and of radius r. For a nonnegative weight function w on M. d , by L p {R d , w), 1 < p < oo, we 
denote the usual Lebesgue spaces related to the measure dw(x) = w(x)dx (in the sequel 
we will often abuse slightly the notation and use the same symbol w to denote the measure 
induced by a density w). Writing X < Y indicates that X < CY with a positive constant 
C independent of significant quantities. We shall write X ~ Y when X < Y and Y < X. 



2. The general setting 

Similarly to numerous frameworks discussed in the literature (see, for instance, the 
setting of deformed Fock spaces discussed by Lust-Piquard [5]) it is reasonable to define, 
at least formally, the Riesz transform lZ k = {Rf, ■ ■ ■ , R%) associated with = — Afe + \\x\\ 2 

as 



(2.1) ld = Sj{C. 



-1/2 



where is a suitable self-adjoint extension in L 2 (M. d , u> fc ) of L k , and are appropriately 
defined first order differential-difference operators. In the present setting we define the 
jth partial derivative 8j related to L k by 

8 j = Tf + x j . 

Taking into account the fact that w k is G-invariant, a short calculation shows that the 
(formal) adjoint of 5j in L 2 (R d ,w k ) is 

To be precise, this means that 

(2-2) (8jf,g)Li(R d ,w k ) = (f,$jg)L2{R d ,w k )i f,9 ^ C™ G (R d ). 

One of the facts which motivate the definition (12. lft is that, as a direct computation shows, 

1 d 

(2.3) L k = -J2(^ + ¥j)- 

3=1 

In the setting of general Dunkl's theory Rosier [9] constructed systems of naturally 
associated multivariable generalized Hermite polynomials and Hermite functions. The 
construction starts from an arbitrarily chosen orthonormal basis {(f n } of V equipped with 
the generalized Fischer inner product 

[P,q}k= (p(T fc )g)(0), p,qEV, 
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such that <f n G V\ n \ and the coefficients of (p n 's are real (more precisely, we should consider 
{(f n } to be a basis of the completion of V). Above, p(T k ) is understood as the differential- 
difference operator derived from p(x) by replacing Xj by Tjf. It is clear that V mi -L k V m2 
for mi 7^ m 2 . For polynomials p G V mi , q G V m2 , rri\,mi G N, the following identity holds 
(see [9], (3.1)], actually we use a slight modification of it) 

(2.4) [p, q] k = c -i 2 ( mi+m2 )/ 2 J ^ exp ( - ^)p(x) exp ( - ^ q ( x ) e -^ 2 w k (x) dx, 

where the constant Ck is the so-called Macdonald-Mehta-Selberg integral, 

Cfc = / exp(—\\x\\ 2 )wk(x) dx 

(note that for any pGP the series defining exp(— A^/4)p terminates). 
The generalized Hermite polynomials H k , n G N d , are then defined by 

H k (x) = 2H v / ^!exp(-A fc /4K(x). 

Note that the definition depends essentially on the choice of the system {tp n } which is 
not specified above. Orthogonality of y? n 's and H2.4[) then show that {H k : n G N d }, is an 
orthogonal system in L 2 (R d ,e-^ 2 w k ), while the system of generalized Hermite functions 

h k n (x) = (^^-) 1/2 exp(-||x|| 2 /2)^(x), x G M. d , nE N d , 

is an orthonormal basis in L 2 (M> d ,Wk), cf. [9, Corollary 3.5 (ii)] (note a slightly different 
normalization: adjusting the factor \fn\ in the definition of H k and the canonical choice 
of the basis (f n (x) = (n\)~ 1 l 2 x n results in obtaining the classical Hermite polynomials for 
k = 0, whereas the coefficient in the definition of h^ n normalizes the system in L 2 (M. d , Wk))- 
Moreover, are eigenf unctions of L^, 

L k h k = (2\n\ + 2 1 + d)h k , 

where 7 = F° r k = and (p n (x) = (n!)~ 1 / 2 x n , H®(x) = n*=i H nj (xj), where 

H n p n j £ N, denote the classical Hermite polynomials of degree rij, cf. |9, Example 3.3 
(1)]; similarly, /i° are the usual multi-dimensional Hermite functions, cf., for instance, [15] . 
Let (■, -)k be the canonical inner product in L 2 (R d , w k ), 

(f,9)k= f(x)g(x)w k (x)dx, 
Jw. d 

and || • \\k be the norm induced by (•, -)fc- The operator 

C k f= X)(2|n| + 2 7 + d)</,/l*> fc /l*, 

neN d 

defined on the domain 

Dom(A) = {/ G L 2 (R d ,w k ) : ^ l( 2 M + 2 7 + d)(f, h k n ) k \ 2 < 00}, 

neN d 
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is a self-adjoint extension of L k considered on C^° G (M. d ) as the natural domain (the inclusion 
C^° G {R d ) C Dom(£ i .) may be easily verified). The spectrum of C k is the discrete set 
{2m + 27 + d : m £ N}, and the spectral decomposition of C k is 

00 

C k f = ( 2m + 2 7 + d)Kf, f e Dom(£ fe ) , 

m=0 

where the spectral projections are 

|n|=m 

Since the spectrum of is separated from zero, 

00 

c k 1/2 f = J2(^ + ^ + d)' 1/2 Kf 

m=0 

is a bounded operator on L 2 (M. d , w k ). It would be tempting to furnish the rigorous defini- 
tion of Rj as an L 2 -bounded operator by writing 

(2.5) R)f = (2M + 27 + dy^if, h k n ) k 6jh k n , 

nGN d 

for / with the expansion / = J2n£N d (fi ^n)*^n- This * s indeed possible when two addi- 
tional hypotheses are satisfied: 

(1) the system {Sjh 1 ^ : n £ N d }, after excluding null functions, is an orthogonal system 
in L 2 (R d ,w k ); 

(2) ||5XIU<(H + 1) 1/2 , neN d . 

It occurs that the required properties are consequences of hypotheses that can be imposed 
on {ip n : n £ N d }. For the purpose of the proposition that follows, for j = 1, . . . , d, we 
denote by N d the set of those n £ ~N d for which T k ip n does not vanish identically; also, by 
|| • \\[ t ] k we denote the norm in V induced by [•, -} k . 

Proposition 2.1. Let {ip n : n £ N d } be an orthonormal basis in V such that ip n £ V\ n \ 
and the coefficients of ip n 's are real. Assume, in addition, that for every j = 1, . . . , d, 

(i) the system {T k (p n : n £ N d } is orthogonal in (V, [-, •]&); 

(ii) \\T*<p n \\ [t]k < (|n| + n G NJ. 

Then, for every j = 1, . . . , d, R k defined by (12. 5p is bounded on L 2 (R d ,w k ). 

Proof. Assumptions which we impose on {Tj(p n } imply the corresponding properties of 
{Sjh^} listed above. Indeed, since commutes with A k , thus also with exp(— A fc /4), and 
Tj(fg) = {T]ff)g + (T k g)f provided / is G-invariant, we have 

Sjhi = (c k 2^n\y l/2 e^ 2 / 2 T'H k n = {2^ / c k ) 1/2 e^ ' 2 exp(-A fc /4)(7}V»)- 
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Therefore (1) and (2) follow by applying (12.41) to the polynomials T^ip n . In consequence 
the conclusion follows by Bessel's inequality applied to the orthogonal system {T^{p n : n G 
Nj}. □ 

A comment is probably in order, concerning a possibility of constructing an orthonormal 
basis with properties as in the first sentence of Proposition 12. II and satisfying, in addition, 
at least Condition (i) of that proposition. As indicated in 0, Section 3], since V mi -Lfc V m2 
for mi 7^ 777,2, 777,1, 777,2 G N, the desired basis {(p n } can be constructed by Gram-Schmidt 
orthogonalization within every V m , m G N, from an arbitrarily ordered real-coefficient 
basis of V m . To guarantee fulfillment of Condition (i) in Proposition 12.11 again it is 
sufficient to act with every m > 1 separately; this is because T^V m C V m -i- Summing 
up, given m > 1, an orthonormal real-coefficient basis {(p n : \n\ = m} of "P m should be 
chosen in such a way that for every j = 1, . . .,d, the system {T*fip n : \n\ = m}, after 
excluding null functions, is orthogonal in V m -\- Note that this means that for every 
j — 1, . . . , d, a number of functions T^(p n , \n\ = m, must vanish. 

We mention at this point that in the case of the group Z, d discussed in detail in the 
subsequent sections, the choice of appropriately normalized monomials x n , n G N d , leads 
to an orthonormal basis in V satisfying the assumptions of Proposition 12.11 This may 
be verified directly (see a comment in Section [3]), but taking into account the explicit 
form of the generalized Hermite functions that appear, it may be checked as well that the 
properties (1) and (2) preceding Proposition 12.11 are also satisfied (see Section 0] below) . 
Whether or not such a choice of <^ n 's may be done for an arbitrary reflection group remains 
an open question. 

Another remark that seems to be worth of mentioning is that it would be tempting 
to justify the boundedness on L 2 (IR d , Wk) of Rj given by (12.11) . for an arbitrary reflection 
group G, without requirements (1) and (2), but instead with using an argument based 
merely on (12.11) . H2. 21) and the self-adjointness of With the notation Rj = b*L^^ 2 the 
argument would proceed as follows: 

Ki&<l|j#/l& + n^/i& 

= (5;5 1 £~ k l/2 f,C~ k l/2 f) k + {5^C- k l/2 fXl l/2 f) k 

<((E(^+^)A- 1/2 /),^ 1/2 /) fc 

i=i 

= 2(L k £ k 1/2 f,£ k 1/2 f) k = 2\\f\\h- 

Note, however, that the use of (I2.2p is limited to a special class of functions that would 
require / from a dense subspace of L 2 (IR d , to fc ) (actually from a subspace of Dom(£ fc ^ 2 )) 
to satisfy: C k 1 ^ 2 f ,5jC k 1 ^ 2 f G C^° G (M. d ). It seems to be far from obvious to determine 
whether such a dense subspace always exists. 
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3. Preliminaries to the Z d group case 
Consider the finite reflection group generated by aj, j = 1, ... , d, 

<Tj {x\ , . . . , Xj , . . . , Xd) {x\ , . . . , Xj , . . . , X(i) , 

and isomorphic to = {0, l} d . The reflection Oj is in the hyperplane orthogonal to e 
the jth coordinate vector in R d . Thus R = {±\^2ej : j = 1, . . . , d}, R + = {V%ej '■ j = 
1, . . . , d}, and for a nonnegative multiplicity function k: R — > [0, oo) which is Zg-invariant 
only values of k on R + are essential. Hence we may think k = (ai + 1/2, . . . , + 1/2), 
aj > —1/2. We write aj + 1/2 in place of seemingly more appropriate aj since, for the 
sake of clarity, it is convenient for us to stick to the notation used in [8]. 

In what follows the symbols T°, Sj, A a , w a , L a , h®, [■, -] a an d so on, denote the objects 
introduced in Section [2] and related to the present setting. Thus the Dunkl differential- 
difference operators T®, j = 1, . . . , d, are now given by 

I?f(x) = djf(x) + (a, + 1/2)^^^, / G C\R% 



and the explicit form of the Dunkl Laplacian is 

A *f^ = E [dx-W + —^-^) - («, + V2) ~, ) ■ 

Note that A a , when restricted to the subspace 

C\R d ) R = {fe C\R d ) :Vj 1 d. f(x) = f(a jX )}, 

coincides with the multi-dimensional Bessel differential operator ^f =1 (6 | | + —0-dj), and 
L a = —A a + || x || 2 reduces to 

h ll9 2a;,- + 19 

-A+||x|| -\ 



J=l J J 

the operator investigated in [8] (to be precise, both operators are assumed to act on 
functions defined on R\ = (0, oo) d ). 

The corresponding weight w a has the form 

d 

w a (x) = ]J\x j \ 2a i +1 ~ Yl \((3,x) a \ 2k ^\ xeR d ; 

here (■, •)„ denotes the inner product in L 2 (M. d , w a ). In dimension one (see [HI Example 
3.3 (2)]) for the reflection group Z 2 and the multiplicity parameter a + 1/2, a > —1/2, 
the polynomial basis {y? n } (orthonormal with respect to [•, -] a ) is determined uniquely (up 
to sign changes) by suitable normalization of the monomials x n , n 6 N. Let a n)Q , be the 
sequence determined by the following recurrence relation: for n > 1, a n Q , = na„_ ljQ , if 
n is even and a n ^ a = (n + 2a + l)a n _i ia if n is odd, ao, a = 1. It is easily checked that 
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Vn( x ) = a "-J x n is an orthonormal basis. Moreover, one obtains as the corresponding 
generalized Hermite polynomials the genuine generalized Hermite polynomials Hn +1 ^ 2 on 
K, as defined and studied in pQ, see also [IT] (this is the only inconsistency in our notation: 
writing H% would collide with the already accepted notation). The system is orthogonal 
in L 2 (R, \x\ 2a+1 e~ x2 dx), and 

(3.1) H«: 1/2 (x) = (-l) n 2 2n n\L a n (x 2 ), 

H%U 2 {x) = {-lT2^n\xL a n + \x 2 )- 

here for a > — 1 and n G N, I™ denotes the Laguerre polynomial of degree n and order a, 
see [U p. 76]. Note that when a = —1/2, H® is indeed the classical Hermite polynomial, 
see [U p. 81]. The <i-dimensional setting emerges through a tensorization one-dimensional 
settings. Thus, given a = (cti, . . . , a^) G [—1/2, oo) d , and taking <fin( x ) = a n}d 2 x n , n G N d , 
a n,a = rii=i a ni,ai, the corresponding generalized Hermite functions are given by 

K( x ) = h%(xi) -...-h^xa), x= (x 1 ,...,x d ) eR d , n = {n u ...,n d )eN d ; 

here are the one-dimensional functions 

with the normalization constants c nuCei , where c^m,^ = (2 Ani T{jii + l)r(n« + «j + l)) -1 ^ 2 
and c 2ni+ i, a . = (2 4ni+2 r(n* + l)T{ ni + an + 2))" 1 / 2 . In view of d3U), one has 

lb 2ni+l\ x i) ~ u 2n 4 +l,a, ; e X l U n . {X i ), 

where 

/ r(m + i) \ 1/2 ( r(ni + i) \ 1/2 

d2ni > ai = { - 1)n \r(n i + a i + l)) ' d2ni+1 ' ai = { ~ ir \r(n i + a i + 2)J ' 

Let us comment at this point that {tp^ : n G N d } satisfies assumptions (i), (ii) of 
Proposition 12. 1[ Indeed, T"<£>" = b n a x n ~ ej , where & n>Q , equals a^i 2 multiplied either by 
rij or rij + 2a j + 1, depending on whether rij is even or odd, respectively. Hence (i) holds. 
On the other hand, ||T^Vnll [■,■]<* = ^n,a a n 2 ej , a = C(M 1//2 ), which implies (ii). 

The system {h% : n G N d } is an orthonormal basis in L 2 (M. d , w a ) consisting of eigen- 
functions of L a , 

L a h a n = (2\n\ + 2\a\ + 2d)h«, 

where by \a\ we denote \a\ = a\ + . . . + a d (thus |a| may be negative). For a = 
( — 1/2,...,— 1/2) we obtain the usual Hermite functions. Recall that, considered on ap- 
propriate domain, the operator L a is positive and symmetric in L 2 (M. d ,w a ). The spectral 
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decomposition of C a is 

oo 

C a f = ^{2m + 2\a\+2d)V^f, f G Dom(£ a ), 

m=0 

with the spectral projections 

\n\=m 

The semigroup T" = exp(— tC a ), t > 0, generated by C a is a strongly continuous 
semigroup of contractions on L 2 (R d ,w a ). By the spectral theorem, 



(3.2) T t a f = e- t{2m+2M+2d) V«f, feL 

The integral representation of T" on L 2 (R d ,w a ) is 

(3.3) T?f(x) = [ G?(x,y)f(y) dw a (y), x eR d 



where the heat kernel is given by 



oo 



(3.4) G?(x, y) = J2 e-'( 2m+2 ' Q l +M ) £ K{x)K{y). 

m=0 \n\=m 

In dimension one, for a > — 1/2 it is known (see, for instance, [9j Theorem 3.12] and 
P p. 523]) that 



G °( I ^) = (2ii^iF exp (-5 C0th ( 2t »( j:2+ i' 2 ) 



(xy) a (xy) a+1 



with being the modified Bessel function of the first kind and order u, 

h{z) = ^ r(fc + i)r(fc + ^ + i) ' 1 argz| < ^ 

Jy, as a function on IR + , is real, positive and smooth for any v > —1, see jH Chapter 5]. 
Notice that the ratio I u (z)/z u occurring in the formula for Gf(x,y) represents the entire 
function of z. 

Therefore, in d dimensions, 

G?(x,y)= G t' e (^vl 

ee{0,l} d 



where 
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By means of the standard asymptotics for I v (see [H (5.16.4), (5.16.5)]) it follows that the 
integral in (13.31) converges for any / G L p (R d , w a ), 1 < p < oo, and therefore it may serve 
definition of T"f on these spaces. 
Note that for e Q = (0, . . . , 0) the component kernel G"' £o (x,y) coincides, up to the 
factor 2~ d , with the heat kernel associated with Laguerre function expansions studied in 
[8]. Moreover, by Soni's inequality [12J, 

I„+i(z) < I v (z), z>0, v>-~, 

it follows immediately that for a G [— 1/2, oo) d 

< G?(x,y) < G?' £ °(x,y), t > 0, x,y G R d . 

Then, as a consequence of [El Theorem 2.1], we get an interesting result on the maximal 
operator Tf*f = sup t>0 \T"f\, which may be of independent interest (for the definition of 
the weight classes A", see Section 0] below) . 

Theorem 3.1. Let a G [-l/2,oo) d . Then T? is bounded on L p (R d ,Wdw a ), W G A%, 
1 < p < oo, and from L x (]R d , Wdw a ) to L l '°°(M. d , Wdw a ), W G Af. 

The proof of Theorem 13.11 reduces to an application of the abovementioned result from 
[8]. To see that this is possible, one has to take into account that G^' £o (x,y) is even with 
respect to each coordinate of (x, y) and that w a (x) has the analogous property. Also, it is 
necessary to know that, given a G [—1/2, oo) d and 1 < p < oo, the condition W G A® is 
equivalent to the condition G A^ ,+ for all £ G { — 1, l} d , where W% denotes the weight 
on given by W^(x) = W(£x), and A^ ,+ stands for the Muckenhoupt class of weights 
associated with the homogeneous space || • ||), being the restriction of w a to 

M^; see (1.3), (1.4) and related comments in [8]. We leave further details to an interested 
reader. 

Although in the present section we focused on the 1> d case, it should be pointed out that 
there is a general background for the results considered here for an arbitrary reflection 
group, see [10] for a comprehensive account. In particular, the heat (or Mehler) kernel 
(13. 4p has always a closed form involving the so-called Dunkl kernel, and is always strictly 
positive. This implies that the corresponding semigroup is contractive on L°°(M d ,w fc ), and 
as its generator is self-adjoint and positive in L 2 (M. d , Wk), the semigroup is also contractive 
on the latter space. Hence, by duality and interpolation, it is in fact contractive on all 
L p (R d ,w k ), 1 < p < oo. 

4. Z d RlESZ TRANSFORMS 

First of all, we shall see how Sj's act on h°. It is sufficient to consider the one-dimensional 
situation and then distinguish between the even and odd cases. Recall that 5j = T? + Xj] 



RIESZ TRANSFORMS FOR THE DUNKL HARMONIC OSCILLATOR 11 

in the one-dimensional case we simply write S in place of 5\. For n G N and a > —1/2, 
combining the fact that h% n is an even function with the identity 

(4-1) Tx L - = - L ^ n " X ' 

see [U (4.18.6)], one easily obtains 



Shi, = VAnh^. 

Similarly, using the fact that fof n+ i i s an °dd function, (14. ip and the identities 
(4-2) - ylfcliy) + (a + l)L a n +_\{y) = nL«(y), n > 1, 

which can be deduced from ( l4~Tj) and [U (4.18.2), (4.18.7)], and 
(4-3) U£\ - L% +1 = -L« n>l 

(this is [H (4.18.5)] after correcting a misprint in one of the superscripts), one gets 



5h% n+1 = V±n + 4a + 4h% n . 
Summarizing, in d dimensions, for n G N d and a G [—1/2, oo) d we have 
(4.4) 5 3 h a n = m( nj ,a 3 )h a n _ ej , 

where 

v/2ny, if raj is even, 



m(rij, otj) 



^2rij + 4a j + 2, if rij is odd; 



here, and also later on, we use the convention that h n _ ej = if =0. Note that for 
a = (-1/2, . . . , -1/2) this is consistent with [13, (3.2)]. Notice also that {5 3 h* : n, > 1} 
is an orthogonal system in L 2 (M d ,w a ) and ||<5j/i"|| = 0(|n| 1//2 ), hence the hypotheses (1) 
and (2) in Section [2] are satisfied. 

We now provide a rigorous definition of the Riesz transforms Rj, j = l,...,d, on 
L 2 (R d ,w a ) by setting 

(4-5) R«f = £ m™}^ n M > K)a K_ ep 

for / with the expansion / = J2 n en d (fi K^aKl- Obviously, is bounded on L 2 {R d ,w a ). 
For j — 1, . . . , d, we define the Riesz kernels Rj(x, y) as 

1 f°° 

(4.6) R?(x,y) = -= b 3 . x G a t (x,y)t- x l 2 dt. 

V 71 " Jo 

The following result shows that the kernel R"{x, y) is associated, in the Calderon-Zygmund 
theory sense, with the operator R" defined in (14.51) . 
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Proposition 4.1. Assume that a = (ax, . . . , a^) is a multi-index such that cti > — 1/2 
and Rf is defined by (14.51) . Let f,g£ C™ zd (M. d ) have disjoint supports. Then 



(4.7) (R«f,g) a = / / R°(x,y)f(y)g(x)dw a (y)dw a (dx). 

JR d JR d 

Proof. Let / = E„ e ^ < h n, 9 = E„eN^ KK- Th en 



R-f = E , ===== a£h 
^ \/2h +2 a +2d 



a h a 

n n—ej 



(convergence of the above series is in L ) and, by Parseval's identity. 

m(nj,OLj) 



( 4 -8) g)a = E , 9 / :^ 

„ti V / 2|^l+2|«| + 2 C / 

To finish the proof it is now sufficient to observe that the right sides of (14. 7p and (j4.8p 
coincide. This follows by noting that the proof of (14. 9 j) in Theorem 14.21 below contains a 
proof of the slightly stronger estimate (see Section [5]) 

1 

w a (B(x, \\y - x\\)) 
and thus the assumption made on the supports of / and g implies 



|^G t Q (a;,y)k- 1/2 rft< 



\S jiX Gf(x,y)\t 1/2 dt\g(x)f{y)\dw a {y)dw a (x) < 



oo. 



At this point the desired conclusion easily follows by repeating arguments from the proof 
of [131 Proposition 3.2]. □ 

The next theorem is the heart of the matter: it says that the kernels of Riesz transforms 
associated with the Dunkl harmonic oscillator satisfy standard estimates in the sense of 
the homogeneous space {R d , w a , \\ ■ ||). The corresponding proof is rather long and tricky, 
hence we decided to locate it in a separate unit (Section [5] below). 

Theorem 4.2. Let a G [—1/2, oo) d . Then the kernels Rf(x, y), j — 1, . . . , d, satisfy 

(4-9) \R?(x,y)\< m( \ rprr j x^y, 

3 w a (B(x, \\y-x\\)) 

\\V x , y R^(x,y)\\ < — — - \ —, x^y. 

J \\x - y\\ w a (B(x, \\y - x\\)) 

It is well known that the classical weighted Calderon-Zygmund theory works, with 
appropriate adjustments, when the underlying space is of homogeneous type. Thus we 
shall use properly adjusted facts from the classic Calderon-Zygmund theory (presented, 
for instance, in [3]) in the setting of the homogeneous space (R d , w Q , \\ ■ ||) without further 
comments. 
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For 1 < p < oo, we denote by A" 
related to the space (M d ,w; Q , || ■ ||). 
functions W G L} r 



sup 

BeB 



1 



= Ap (IR d , w a ) the Muckenhoupt class of A p weights 
More precisely, A^ is the class of all nonnegative 

and 



such that W -1 /(p- 1 ) G L\ oc (R d , w. 



w Q (B) 
when 1 < p < oo, or 



W(x) dw a (x) 



B 



w a (B) 



dw a {x) 



B 



p-1 



< OO, 



If 1 

sup — — / W(x) dw a (x) ess sup— — — 
BeB w a {B) J B X&B W[x) 



< oo, 



if p = 1; here £> is the class of all Euclidean balls in M. d . It is easy to check that (in 
dimension one) the power weight W r (x) = \x\ r belongs to A" 1 < p < oo, if and only if 
-(2a + 2) < r < (2a + 2)(p - 1), and W r G A? if and only if -(2a + 2) < r < 0. 

Theorem 4.3. Assume that a = (a±, . . . ,ad) is a multi-index such that ai > —1/2, 
i — 1 d. T/ien the Riesz operators R" , j = 1, . . . , d, defined on L 2 (M d , w a ) by i \4-5§ , 



are Calderon-Zygmund operators associated with the kernels R^(x,y) defined by In 
consequence, R® extend uniquely to bounded linear operators on L p (M d ,Wdw a ) , 1 < p < 
oo, W G A^, and to bounded linear operators from L 1 (R d , Wdw a ) to L 1,00 (IR d , Wdw a ), 
W G A\. Denoting these extensions by the same symbols R?, for f G L p (R d ,Wdw a ), 
1 < p < oo, W G A p , we have 



(4.10) 



m(nj, aj) 



n\ 



2\a\ +2d 



{f,K) a . 



Proof. The first statement is just a combination of Theorem 14.21 and Proposition 14. 11 The 
second is a consequence of the general theory, see [3]. To verify (14.101) fix j = 1, . . . , d, 
and choose f k G L p (R d , Wdw a ) n L 2 (R d ,w a ) such that f k ^fm L p (R d ,Wdw a ). Then, 
for any n G N d , (f k ,h%) a - (/,/£)„ and (Rp k ,h%_ e .) a - (R«f, h%_ e .) a , where Rff is, 
by the very definition, the limit of Rffk in L p (R d , Wdw a ). The claim follows. □ 



We finish this section with an application of Theorem 14.31 connected with a priori 
L p -bounds. This kind of estimates is a useful tool in the theory of partial differential 
equations. 



Proposition 4.4. Let a G [—1/2, oo) d , 1 < p < oo, and i,j G {l,...,d}. Then 

,w a ) ^ W-^aJ \\LP(R d ,w a ), 



\$i3jf\\LP(R d ,w a ) < \\L a f\\LP(K<i,w a ), f ^ C™"-"" 



Proof. Let Rf be the adjoint of Rf taken in L? 



Then 



m(rij + 1, aj) 



neN d 



+ 2|ct| +2d + 2 



feL* 



w a ). 
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By Theorem 14. 31 and duality, R® extends to a bounded operator in L p (M. d , w a ), 1 < p < oo. 
Therefore, to finish the proof, in suffices to show that 

5*5, f = R?R«L a f 

for / G C^°(WL d ). However, since such functions can be suitably approximated by linear 
combinations of h^, it is enough to verify the relation for / = h°. 
Using the identity 

StK = m{m + l,a i )h^ e ., neN d , 

which, taking into account that 5* = —5i + 2xj, can be obtained with the aid of (14.41) . 
f l42|) and we get 



5*5jh® = m(rij, aj) 5*h"_ e . = m(nj, m((n — e 3 -)j + 1, «;) h"_ 
On the other hand, 

R?R«L a K = {2\n\+2\a\+2d)R?( m ( n ^ a ^ u 



" ! ' ' ! ' V v/2|n|+2|a| + 2|d| 3 

m{rij,0ij) m((n — ej)i + l,ai) 



= (2\n\ + 2\a\ + 2d)- / ,_. , 

V^2|n| + 2\a\ + 2\d\ y/2\n - e s \ + 2\a\ + 2d + 2 3+1 

= m[n h ay) m({n - e s )i + 1, on) K_ e . +e .. 
The conclusion follows. □ 



5. Kernel estimates for the Zfj Riesz operators 

This section is devoted to the proof of the standard estimates stated in Theorem 14. 21 
We begin with a suitable decomposition of the kernels under consideration: 

£G{0,l} d 

where the component kernels are given by 

1 f°° 

(5.1) Rf\x,y) = -= SjGFfay^-Wdt; 

v n Jo 

here and also later on 5j always refers to the x variable. Clearly, it is sufficient to prove 
the estimates in question for each R" ,e (x,y), e G {0, l} d , separately. Thus, from now on, 
we fix e G {0, l} d . 

Next, we claim that it is enough to show the relevant bounds only for x,y G IR+ (here, 
in fact, we could replace by any other Weyl chamber of the considered group of 
reflections). Observe that G^ ,£ (x,y), as a function of either x or y, is even with respect 
to the jth axis if Sj = and is odd with respect to the jth axis if Sj = 1. The derivative 
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8j, when taken with respect to x and restricted to even functions in the jth coordinate, 
has the form 

(5.2) S=£ + ** 
whereas the restriction of 8j to odd functions in the jth coordinate is 
5.3 5°= — + Xj + ^ . 

Thus 5jG"' £ (x,y), as a function of x, is either even or odd with respect to the jth axis, 
depending on whether Ej = 1 or Ej = 0, respectively. Moreover, as a function of y, it is 
either even or odd with respect to the jth axis, depending on whether Ej — or Ej — 1, 
correspondingly. It follows that 

\R? e (x,y)\ = \R? £ ( V x,Zy)\, v ,t G 

Similarly, we have 

\\V x , y Rf £ (x,y)\\ = \\V x>y R^(r]x^y)l 77, £ G 
This, together with the symmetry of w a , 

w a {x) = w a (£x), £G{-l,l} d , 

proves the claim. 

Finally, we also observe that, after restricting x and y to R+, with no loss of general- 
ity we may replace the Euclidean balls B(x, \\y — x\\) in the relevant estimates by their 
intersections with Rl, denoted further by B + (x, \\y — x\\). 

Summing up, we are reduced to proving the following. 

Lemma 5.1. Assume that a G [-l/2,oo) d , e G {0, l} d and j G {l,...,d}. Then 
\R? e (x,y)\< 1 — , x,yeR d + , x^y, 

\\V x , y Rf e (x,y)\\ < 1 1 — , x,yeR d + , x^y, 

\\x-y\\w+{B+(x,\\y -x\\)) 

with w+ being the restriction of w a to R+ . 

The proof of Lemma 15.11 relies on a generalization of the technique applied recently 
by the authors in the Laguerre function setting [8]. Thus the main tool in performing 
necessary kernel estimates will be Schlafli's integral representation of Poisson's type for 
the modified Bessel function, 



I v (z) = z v J exp(—zs) U u (ds), 



z > 0, v > --, 
- 2' 



i n 
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where the measure ITj, is determined by the density 

TlJds) = ±= y — , se -1,1 

vy ' ^2"T{v+ 1/2)' V 

when v > —1/2, and in the limiting case of v = —1/2, 

n_i/2 = -7= (?7-i + ?7i) , 
V In 



where r]_i and rji denote point masses at —1 and 1, respectively. 

Consequently, for a G [—1/2, oo) d , the kernel G^' £ (x,y) may be expressed as 

(xy) £ f ( 1 



J exp (-icoth(2 t )(|Nr + IMD-Ei^l)) W*l 



2 rf (sinh(2t)) rf +l«l+l £ l i hlil]d 2 /v " y " y ^sinh(2t) 

where n a denotes the product measure ®f =1 n a . . Then the change of variable t 
(5-4) t = t( C ) = ilogi±|, Ce(0,i), 

leads to the following useful symmetric formula: 

1/1 — ( 2 \d+H + \e\ f ( 1 £ \ 

GT £ (£,2/) = 2 ( J (^) e y 11]d exp {~^+( x ^ s )~-^-( x ^ s )) n a+e(^) ; 

with 

q±{x,y,s) = \\x\\ 2 + \\y\\ 2 ±2^2xiyiSi. 

i=i 

For the sake of brevity we shall often write shortly q + or q_ omitting the arguments. 

We pass to the Riesz kernels. By the change of variable (15.41) the component kernels 
(15. ip can be expressed as 

r r 1 

yCt,E I 



(5.5) Rf E (x,y)= U a+£ (ds) (3d, a +e(() Sjipfa, y, s) d(, 

J[-i,i] d Jo 

where 

r tt\ ^ a-c 2 V +|q| 1 (, i+cv 1/2 

(this is, up to the constant factor 2~ d , /3] a (C) defined in [SJ, (5.4)]), and 
^c{x,y,s) = (xy fexp ( - — q + (x,y,s) - -q-(x,y, s)J; 

passing with 5j under the integral sign and then changing the order of integration that 
was implicitly performed above will be justified in a moment. 



' This useful change of variable was invented by Stefano Meda, in the context of the Ornstein-Uhlenbeck 
operator. 
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To compute 6jip^(x, y, s), the expression appearing in (15. 5p . observe that the derivative 
5j may be replaced either by 5| or 5°, depending on whether Sj — or Sj — 1, respectively; 
see ( 15.21) . (15.31) and the accompanying comments. Then one easily arrives at 



(5.6) 6jil%(x,y,s) 



(xyY^Xj - —(.r, ■ y jSj ) 
+ X{e 1 =i}(2a j + 2)y ] (xyY 



c 



(xj yjSj) 



exp \ - — ^ 



We now come back to explaining the possibility of exchanging Sj with the integral sign. 
Observe that integrating the expression 



Ii a+£ (ds) d x 



[-1,1]* 



1 t 
(xy) £ exp[ - —q + {x,y,s) - -q-{x,y,s) 



with respect to Xj in an interval [0,v], < v < oo, leads to an absolutely convergent 
double integral. Thus we may apply Fubini's theorem and then differentiate the resulting 
equality in v. This gives the desired conclusion. 

The application of Fubini's theorem that was necessary to get (15.51) is also justified 
since, in fact, the proof (to be given below) of the first estimate in Lemma 15.11 contains 
the proof of 



[-1,1]" 



U a+£ (ds) / Pd, a +e(() \6jip%(x,y,s)\d( < oo. 
Jo 



For proving Lemma 15.11 we will need several simple technical results obtained in [H 
Section 5]. We gather them in the lemma below for the sake of reader's convenience. 

Lemma 5.2 ([8j Corollary 5.2, Lemma 5.5]). Assume that a G [— 1/2, oo) d . Let b > 
and c > be fixed. Then, for any j — 1, . . . , d, we have 



< c b l 2 



(a) ({xj + yjSjl + \yj + a^l) 6 exp ( - c^q + ( 

(b) ({xj -yjSjl + \yj -XjSj\) b exp c(q„(x,y, s)) <C b/2 , 
c^q+(x, y, s) - c(q-(x, y, s)) < C &/2 , 



[Xj) exp 



(d) 



/MOC 



-6-1/2 



exp 



c^q+(x,y,s)^ dC < (q+(x,y,s)) d H b , 



uniformly in x,y G , s G [—1, l] d and, except of (d), in ( G (0, 1). 

We will also need the following generalization of [8, Proposition 5.9]. 

Lemma 5.3. Assume that a G [—1/2, oo) d and let 7, k G [0, oo) d be fixed. If a complex- 
valued kernel K(x,y) defined on Ml x Ml\{(x,?/) : x — y} satisfies 



\K(x,y)\ <{x + y) 



27 



n 



•i,i] c 



a+7+ 



K {ds) (q + (x,y,s)) 



-d-M-M 
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then also 

\K(x,y)\ < 1 jj^r, x^y. 

w+(B+(x, \\y-x\\)) 



Similarly, the estimate 



\\V x , y K(x,y)\\ < (x + yf [ U a+J+K (ds) (q + (x, y, s ))-H«l-H-i/» 

J [-1,1]* 

implies 

\\V Xt yK(x,y)\\ < j. — 17 1 rrrr, x^y. 

\\x-y\\ w+{B+{x, \\y-x\\)) 

Proof. When n = (0, . . . , 0), the proof is a direct modification of the arguments used in 
[8] for showing Proposition 5.9 there. The difference is that now [H Lemma 5.8] must 
be applied with different parameters, but otherwise the reasoning is essentially the same. 
Thus we omit the details. 

The general case of arbitrary k follows easily from the following observation: for any 
fixed u, v > and j = 1, . . . ,d, 

/ U a+y+U£j (ds) (q + (x,y,s))~ v < U a+J {ds) (q+(x,y,s))~ v . 

J[-i,i[ d J[-i,i] d 

Thus it remains to justify the above estimate. 

When atj + jj > — 1/2 the relation between the (one-dimensional) densities, 

n a j+7j+u(s.j) ^ n Qj . +7i (sj), Sj e (-1, i), 

deduced immediately from the exact formula, does the job. In the case ctj+jj = —1/2 it is 
enough to notice that (q + )~ v , as a one-dimensional function of Sj G [—1,1], is decreasing; 
in particular, it attains its maximum at the endpoint Sj = —1. Then the conclusion follows 
by the special form of H-1/2 and the finiteness of the measure n a . +7 . +u . □ 

Finally, notice that for A G [—1/2, oo) d , j G {1, . . . , d} and u > 0, 

(5.7) /?d,A +uej (C)<r u ^,A(C), Ce(o,i). 

This property follows directly from the exact formula for /3d, a - 
Proof of Lemma \5. 11 the growth estimate. By (15.51) and (15. 6p 

\Rj' £ (x,y)\ ^ [ _ U a+£ (ds) j (3 d ^ +£ {C){xy) £ Xj + ^\xj + yjSj] + (\xj - yjs. 



f-i,i] d ^0 



exp(--g + 

n a +e(ds) J (3d, a +e(()yj(xy) £ ^ exp ( - — q + - -q-^jd( 



4C 

+ X{e j =l} j 
■ll+l 2 - 



RIESZ TRANSFORMS FOR THE DUNKL HARMONIC OSCILLATOR 19 

We treat the integrals Ti and X 2 separately. Applying Lemma [5T21 (a), (b) and (c) we get 

1i < (xyf I U a+£ (ds) [ (3 dia+E (C)C 1/2 ex V (-^q+- f <7-W 
J[-i,i] d Jo v °C <=> y 

Next, using Lemma [5.21 (d) and observing that (xy) e < (x + y) 2e gives 



Ii<(x + y) 2£ [ U a+e (ds)(q, 
J\-i,l] d 



>— d— \a\ — \e\ 



Now Lemma [5T3l taken with 7 = 5 and k — (0, . . . , 0), provides the required growth bound 
for Xi . 

To estimate X 2 we assume that Ej = 1 (otherwise X 2 vanishes) and use (15. 7p . specified 
to A = a + £ — ej/2 and u = 1/2, getting 

X 2 < / U a+£ (ds) [% d>a+e . ej/2 (C)(x + y) 2 ^C 1/2 exp(-j-q + -L)dC 
J[-x,i] d Jo v % 4 ' 

Then Lemma 15.21 (d) shows that 

I2 < (x + yf^M [ U a+£ (ds)(q + )- d -^- e ^. 
J[-i,i] d 

Combining the above with Lemma [5.31 applied with 7 = 5 — ej/2 and k = ej/2, produces 
the relevant bound for X 2 . This finishes proving the growth estimate for R®' e (x, y). □ 

Proof of Lemma I5.il the smoothness estimate. By symmetry reasons, without any loss of 
generality we may focus on the case j = 1. Thus the kernel under consideration is 

Ri' £ (x,y)= / U a+£ (ds) / (3 dta+£ (()8 1 il> l ;(x,y,s)d(, 
J[-i,i] d Jo 

with Sii/j^(x, y, s) given by (15.61) . Our task is to estimate || V Ei j / -R"' e (a;, y)\\. However, since 

r r 1 

W^x,yRl \X 



< 



(ds) / (3 dt0i+£ (() V a! , y [(JiV'f(ar,2/,s)] dC, 
i,i] d Jo 

U a+£ (ds) / (3 d>a+£ (C) Hv^^i^^,?/^)]!! d( 

-lA] d Jo 



'[-1,1]" 

(passing with V x ,y under the integral signs is legitimate, the justification being implicitly 
contained in the estimates that follow; see comments at the end of this proof), and since 
only the first variables Xi,yi are distinguished, it is sufficient to show a suitable bound 
for the double integral in (15.81) . but with the gradient restricted only to the first two 
coordinates in x and y. 

To proceed, we need to compute the relevant derivatives. Denoting \l/"' e = 6iip^(x, y, s) 
and abbreviating the exponential entering ip^(x, y, s) to Exp(£, q±) we get 

X{ei=i}Z/i(^) £ ^ ei (xi - ^-{xi + y x si) - ^(xi - yisi)) 



dxi^^ 



20 
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_d_ 

dyi 



X{ e i=i}Xi(x?/) £_ei (xi - ^(xi + yisi) - -(a* - yisi)j - (xy) 



C ' 



£l _ C«l 

2C 2 



+ 2(ai + l)X{«=i}(a^) 



Exp(C,g±) - (^(yi + zisi) + -xisi))#"'' 



dx 2 ^ 



1 ^ 

X{ e2 =i}Z/2(xy) £ ~ e2 (xi - — (xi + yisi) - -(xi - yisi)) + 2{a x + 1) 



X{ei=£2 = 



e— ei— e2 



Exp(C,g±) - (^( x 2 +2/2S 2 ) + ~(x 2 -y 2 s 2 ))^' £ , 



<9j/ 2 



1 ^ 

X{ £2 =i} a; 2(a;i/) e_e2 (x! - — (xi + y x si) - -(x x - j/jsi) ) + 2(«i + 1) 



X{ £l = £2 =i}yiX2(x?/) £ 



Exp(C, g±) - + x 2 s 2 ) + ^(y 2 - x 2 s 2 )) * 



O.c 

c • 



We will estimate separately the four integrals resulting from replacing in (15.81) by one 
of the above derivatives. Denote these integrals by J Xl) J yi , J X2 and Jy 2 , respectively, 
and assume that e\ — e 2 — 1 (otherwise the situation is even easier since some terms do 
not appear). 

Applying (a), (b) and (c) of Lemma [5.21 we get 

J Xl <(x + y) 2£ ~ ei [ n a+ e(ds) [ ^, Q+£ (C)C" 1/ VExp(C,g ± )^C 

J[-l,l] d Jo 

+ (x + y) 2£ [ U a+£ (ds) [ f3 d , a+e (C)C lE MC,q±)dC 
J[-i,i] d Jo 



+ 



/ U a+£ (ds) [ p d , a+£ (()C 1/2 Ml(x,y,s) 
J[-i,i] d Jo 



lExp(C,g±) 



To estimate we use (15.71) and then Lemma [5.21 (d). obtaining 

J*> £ {x + yf^ [ U a+E (ds) I /3 d , a+£ . ei/2 (OC VExp(C, q±) d( 
J[-i,i] d Jo 

<(x + y )^-^ [ n a+£ (ds)( q+ r d -^-^\-y\ 

J[-i,i] d 

The expression can be treated in the same way as X± appearing in the proof of the 



growth estimate above. We get 



Jl < (x + y) 2e [ n Q+e 



{ds)(q + ) 



-d-\a\-\e\-l/2 



RIESZ TRANSFORMS FOR THE DUNKL HARMONIC OSCILLATOR 21 

Finally, estimating is completely analogous to estimating the growth of \R" ,e (x, y)\ 
performed earlier. The result is 

Jl < (x + yf £ [ TL a+e (ds)( g+ r d -^- 1 / 2 
J[-i,i] d 

+ (x + yf^m f Y[ a+£ {ds){q + )- d -^ £ ~ ei/2 \- l/2 . 
J[-i,i] d 

Combining the above estimates of . J~ 2 and with Lemma |5T3| specified to either 
7 = e and k — (0, . . . , 0) or 7 = e — ei/2 and k = ei/2, provides the required smoothness 
bound for J Xl . 

Considering J yil items (a), (b) and (c) of Lemma [5.21 lead to 

Jin < Jl, + Jl + Jl + (x + yf £ - 2ei I H a+£ (ds) f /W(0 Exp(C, q±) d(. 

J[-i,i] d Jo 

Thus it suffices to bound suitably the last term of the sum, which we denote by J}_. 
Making use of ( 15. 7ft and then applying item (d) of Lemma 15.21 gives 

Jy\ < (x + y) 2 ^ [ U a+£ (d S ) f (3 d , a+ e-e 1 (0C 1 ^MC,q±)d( 

J[-i,i] d Jo 

<( x + yfie-e,) f n a+£ (rf S )(g + )-^l Q H £ - ei l- 1 / 2 . 
J[-l,l] d 

Now Lemma 15.31 employed with 7 = e — t\ and k — ex implies the desired bound of J yi . 
Passing to J X21 items (a), (b) and (c) of Lemma I5T21 reveal that 

J* 2 <(x + y) 2£ - e > f U a+£ (ds) f p d>a+£ (()C 1/2 V^M(,<l±)d( 
J[-i,i] d Jo 

+ ( x + y )*-ei-ea f U a+£ (ds)[ (3 d , a+£ (C) ^MC, <Z±) d( 
J[-l,l] d Jo 

+ [ u a+e (ds) [ ^ a+e (0C 1/2 \Sir c (x,y,s)\^2==dC 

J[-i,i] d Jo A/Exp(C,g±) 

^ X2 1 ^ X2 1 ^ X2' 

Treatment of J X2 and J X2 is the same as in case of J Xl and J% . Thus it remains to deal 
with J 2 r Let 7 = e — ei/2 — €2/2. In view of ( 15.71) and Lemma 15721 (d) we have 

Jl <(x + y) 2 ^ [ U a+£ (ds) I Pd^iQC 1 Exp(C, q±) d( 
J[-i,i] d Jo 

<(x + y)^ [ IL a+£ (ds)(q + y d -W-W-V 2 . 
J[-i,i] d 

This, together with Lemma T5.3I taken with k = ei/2 + e2/2 provides the relevant estimate 
of J X2 . 

Eventually, the analysis related to J V2 is analogous to that for J X2 . 
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We now come back to explaining the possibility of exchanging d Xl with the integral over 
[— 1, l] d x (0, 1) (a similar argument works for the remaining derivatives, and considering 
only j = 1 does not affect the generality). To apply Fubini's theorem and then use the 
argument invoked earlier in connection with (15.51) . it is sufficient to show that 

/ I / Il a+£ (ds) \d Xl [5i4> £ c {x,y,s)]\(3 dia+£ {()d()dxi < oc, 
Ju \J[-i,i] d >^o / 

for any < u < v < oo such that yi is not in [u, v] when Xj = yi, 2 < i < d. This estimate, 

however, may be easily obtained with the aid of the bound for J7 ri established above. 

The proof of the smoothness estimate in Lemma 15.11 is complete. □ 
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